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ABSTRACT 

We analyze the ultraviolet spectrum of the cataclysmic variable GK Per at 
maximum light. The flat ultraviolet spectrum in this system requires a truncated inner 
accretion disk and an unusually flat radial temperature profile. This requirement is 
not satisfied by any non-magnetic steady or non-steady disk model. 

We consider a magnetized accretion disk model to explain the ultraviolet spectrum. 
The available data on the white dwarf spin and possible quasi-periodic oscillations 
constrain the magnetic field, B*, and the disk accretion rate, M, to lie along a 
well-defined spin- equilibrium condition (M/10 17 g s" 1 ) ~ 100(6* /10 7 G) 2 . Our self- 
consistent treatment of the magnetic torque on the disk flattens the disk temperature 
distribution outside the disk truncation radius. This modified temperature distribution 
is too steep to explain the UV spectrum for reasonable field strengths. 

X-ray heating is a plausible alternative to magnetic heating in GK Per. We 
estimate that the disk intercepts ~ 5% of the accretion energy in outburst, which 
results in an extra disk luminosity of ~ 5-10 L@. Model spectra of optically thick disks 
are too blue to match observations. The UV spectrum of an optically thick disk with 
an optically thin, X-ray heated corona resembles the observed spectrum. The X-ray 
luminosity observed during the outburst indicates M < 10 18 g s _1 , which is a factor 
of 10 lower than that required to explain the ultraviolet luminosity. Radiation drag on 
material flowing inward along the accretion column lowers the shock temperature and 
reduces the X-ray luminosity. Most of the accretion energy is then radiated at extreme 
ultraviolet wavelengths. 

Subject headings: accretion, accretion disks — binaries: general — magnetic fields — 
stars: individual (GK Per) — stars: magnetic fields — Ultraviolet: stars 
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The old nova GK Persei (Nova Per 1901) is an intriguing cataclysmic variable (CV). It is 
one of a few old novae known to undergo dwarf nova outbursts, albeit with the unusually long 
recurrence time of ~ 400 d (Sabbadin & Bianchini 1983). The system has a coherent 351 s X-ray 
pulsation (Norton et al. 1988, Norton & Watson 1989, Ishida et al. 1992) that has also been 
detected in ultraviolet light curves (Patterson 1991). Longer wavelength optical light curves show 
quasi-periodic oscillations at periods of 360-400 s (Patterson 1991). Although the optical spectrum 
of GK Per is normal compared to other long period CVs, its ultraviolet (UV) spectrum has a flat 
continuum and strong emission lines compared to other dwarf novae. The emission lines remain 
prominent throughout outburst, even though the UV continuum temperature is cool, ~ 10 4 K. 

GK Per has many properties in common with the DQ Her subclass of CVs, which are 
magnetic binaries also known as intermediate polars (Warner 1985). Most of the system's 'unusual' 
properties are then understood if an accretion disk fed by the lobe-filling K0 secondary is truncated 
by a magnetic field well above the white dwarf photosphere. The X-ray pulsations result from 
'hot spots' where the magnetic accretion columns impact a white dwarf with a rotation period of 
351 s. The flat UV spectrum and the long recurrence time for dwarf nova outbursts also support 
this interpretation. For outbursts beginning in the inner region of the accretion disk (e.g. Smak 
1991,1993, Kim et al. 1992), assuming that the effective disk temperature in the inner region is 
constant, the recurrence time scale is roughly the viscous time scale oc R 2 /v where R is the radius 
of the thermally unstable inner region of the disk and v is the kinematic viscosity coefficient (Smak 
1993). Using Smak's (1993) formula for the critical surface density, S cr jt, for the cold branch disk 
and v ~ M/37rS cr jt (e.g. Frank et al. 1992), the recurrence time scale can be estimated as 

where a is the Shakura-Sunyaev viscosity parameter (e.g., Frank et al. 1992), M* is the mass of 
the accreting white dwarf, and M is the mass accretion rate in the inner region of the the cold disk 
(which is different from the mass transfer rate). The recurrence time scale depends very weakly 
on the mass transfer rate (e.g. Cannizzo & Mattei 1992). For a typical quiescent a ~ 0.01, the 
observed long recurrence time scale implies a truncated disk with an inner radius of Rq > 10 10 
cm (Cannizzo & Kenyon 1986, Angelini & Verbunt 1989, Kim et al. 1992). The lack of disk 
material close to the white dwarf also produces a flat UV continuous spectrum by reducing the 
characteristic disk temperature. 

Despite the conceptual simplicity of this interpretation, several observations of GK Per remain 
poorly understood. The X-ray luminosity is a small fraction of the UV luminosity in outburst and 
quiescence as in other DQ Her systems (Patterson 1994). Both the flat UV spectrum and the 
outburst recurrence time scale require inner disk radii large compared to the corotation radius of 
a white dwarf with a 351 s spin period, Ro > 10 10 cm vs R c ~ 8 x 10 9 cm. Accretion should not 
occur in the circumstances. The UV spectrum further requires larger values for both M and Rq 
than allowed by disk instability models for the eruption. 
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In this paper, we consider GK Per as a laboratory to test our understanding of the physics 
of the disk-magnetosphere interaction. Besides disrupting the disk at an inner radius larger than 
the stellar radius, the stellar magnetosphere redistributes angular momentum and energy through 
torques on disk material (Ghosh & Lamb 1979ab). This interaction modifies the disk temperature 
distribution and accretion luminosity in a predictable way that observations can test directly 
(Mauche et al. 1990, Kenyon et al. 1996) We focus on the dissipation in the disk with some 
recent improvements and use emission from the accretion column and the disk as diagnostic tools 
for magnetized accretion disk models. Our goal is to use the observed UV spectrum and various 
X-ray constraints to develop a self-consistent picture for this and other DQ Her system. We begin 
with a discussion of simple disk models in §2, describe and apply a magnetized disk model in §, 
consider X-ray heating of the disk in §4, and conclude with a brief summary and discussion in §5. 



2. The Observed UV Spectrum and Simple Disk Models for GK Per 

The UV spectrum of GK Per at maximum is very unusual compared to other dwarf novae 
in outburst (see also Wu et al. 1989). The dereddened continuum is almost flat in F\ over 
1200-3200 A and closely resembles the continuum of a B8-B9 V star with an effective temperature 
of 10,000-12,000 K (Wu et al. 1989). In contrast, the UV continua of most dwarf novae at 
maximum resemble spectra of much hotter B2-B3 V-III stars having temperatures of ~ 20,000 K 
(la Dous 1991). The UV spectrum of GK Per also contains strong He II, C IV, and N V emission 
lines. These lines are usually strong absorption features in dwarf novae at maximum, although 
some have very weak emission lines (la Dous 1991). The strong emission features in GK Per are 
similar in strength to lines observed in other DQ Her-type stars in the high state (la Dous 1991). 
The He II line in GK Per is, however, exceptionally strong compared to lines observed in most DQ 
Her stars. 

Bianchini & Sabbadin (1983) first showed that standard steady-state disk models could not 
reproduce UV observations of GK Per at minimum. Wu et al. (1989) reached the same conclusions 
for outburst spectra. Both groups showed that disks with an inner hole could roughly match 
observed fluxes. Bianchini & Sabbadin (1983) estimated an inner disk radius of Rq ~ 10 9 cm for 
M ~ 2 x 10 16 g s _1 at minimum; Wu et al. preferred a much larger inner radius, Rq ~ 3 — 5 x 10 10 
cm, in outburst when M « 5 — 20 x 10 19 g s _1 . Wu et al. (1989) remarked that the UV data do 
not constrain Rq and M independently, because any disk model with a maximum temperature of 
~ 10 4 K produces a flat spectrum similar to the observed spectrum. Their set of models has M 
slowly decreasing with increasing Rq to reproduce the observed spectrum. 

To compare model disk spectra with observations, we first verified the Wu et al. (1989) 
measurements using IUE data at maximum light collected from the NSSDC archives and binned 
the F\ data, dereddened by Eb~v = 0.3, in 25 A intervals as in Wu et al. The filled circles in 
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Figure 1 plot our results. The observed fluxes have an intrinisic uncertainty of ±10% as indicated 
by the error bar. The probable error in the reddening, 5Eb~v ~ 0.05, introduces an additional 
error of ±25% at 1300 A and ±10% at 3000 A. We then computed spectra for steady-state disk 
models in which each annulus radiates as a main sequence star (see Cannizzo & Kenyon 1986). 
For the central star, we adopt M* = 0.9 M© (cf. Wu et al. 1989, Watson et al. 1985, Mauche et 
al. 1990, Kim et al. 1992) and R* = 6.22 x 10 8 cm from the white dwarf mass-radius relation 
(Nauenberg 1972, Ritter 1985). The disk has a corotation radius of R c = 7.20 x 10 9 cm for a 
rotational period of = 351 s and an outer radius of R ou t = 1-5 x 10 11 cm appropriate for the 
K0 IV secondary (e.g., Cannizzo & Kenyon 1986). We adopt a disk inclination of i = 60° and a 
distance of 470 pc (e.g., Bianchini & Sabbadin 1983, Crampton et al. 1986). The uncertainty in 
the absolute level of our predicted fluxes is roughly ±50% for reasonable uncertainties in the mass, 
radii, inclination, and distance. The slope of our predicted spectra have higher accuracy, ±20%, 
as long as main sequence stars are a reasonable approximation to the emitted spectrum. We thus 
consider a model successful if it matches the observed fluxes to a factor of 2 and the spectral slope 
to ±25%. 

Figure 1 compares spectra of two steady disk models with the dereddened UV continuum 
data. We computed models for R = 0.9 R c and M = 8.8 x 10 18 g s" 1 (SSI) and R = 5 R c and 
M = 4 x 10 19 g s _1 (SS2). Neither model reproduces the data very well. The first model (SSI) has 
an inner radius near the corotation radius and also has roughly the correct bolometric luminosity, 
L s=s 20 Lq. This model - and any other steady disk model with Ro « R c - is too blue due to the 
high temperatures, ~ 2-3 xlO 4 K, near the corotation radius. The second model (SS2) has the 
correct UV slope but is not luminous enough. More luminous models are hotter than this model 
and fail to match the UV slope. We could match the UV data better with a smaller i for model 
SS2, but observations of binary motion in GK Per preclude i < 45° (Crampton et al. 1986). 

Non-steady disk models also fail to match UV observations of GK Per. Kim et al. (1992) 
developed disk instability models to account for the system's dwarf nova eruptions and matched the 
visual magnitude and recurrence time of typical eruptions. We read the temperature distribution 
for their model 82 from their Figure 8 to compute the 'DP spectrum shown in Figure 1. Our 
calculation reproduces Kim et al.'s spectrum in their Figure 9 to ~ 10% as judged by a visual 
comparison. The spectrum of this model resembles our second steady-state model and also falls 
short of the observations. Although the DI and SS2 models reproduce the observed maximum 
visual magnitude of GK Per, both have i?o 3> R c - This result is contrary to our expectations in a 
magnetically truncated disk model. We expect Rq < R c , which allows disk material to match the 
rotational velocity of the central star and accrete along magnetic field lines. 

The failure of standard disk instability and steady-state disk models suggests that another 
mechanism heats the disk. To produce a flat UV spectrum in the disk, we require a flatter radial 
temperature gradient than in a standard steady-state disk model (e.g., T oc R~ x , with x « 
0.5-0.6). Non-steady disks often have such flat temperature gradients in quiescence, when the disk 
mass grows with time (see Kim et al. 1992). Most non-steady disk models closely approximate 
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steady disks with x « 0.75 during outbursts (see Kim et al. 1992 and references therein). The 
UV spectra of currently available non-steady models are also not luminous enough (Figure 1). 
External heating mechanisms, however, can produce flat temperature gradients. Ghosh & Lamb 
(1979ab) first showed that the magnetic torque at the inner edge of a magnetically truncated disk 
acts to flatten the disk temperature gradient close to the corotation radius. Mauche et al. (1990) 
applied this concept to GK Per and reproduced some aspects of the observations. X-ray heating 
is another external mechanism that can flatten the disk temperature distribution. Although disk 
irradiation from an accretion hot spot or boundary layer is unimportant in most CVs (Wade 
1988; Smak 1989), the large disk in GK Per can intercept and reradiate a significant fraction of 
the luminosity emitted from the central white dwarf. The temperature distribution of a concave, 
illuminated disk can approach T oc R~ l l 2 in some disks (e.g., Kenyon & Hartmann 1987). We will 
first consider magnetic heating in §3 and then develop a simple X-ray heating model in §4. 



3. Disk-Magnetosphere Interaction and Magnetic Heating 



3.1. Magnetic Truncation Radius 



To estimate the magnitude of disk heating in the GK Per disk, we need to derive the disk 
truncation radius, Rq, and the magnetic torque on the disk, N. To compute these quantities, we 
follow Yi (1995) and adopt a geometrically thin, optically thick accretion disk in which the viscous 
angular momentum transport is parametrized by the conventional a- viscosity (e.g., Frank et al. 
1992; see also Wang 1987; Campbell 1992; Kenyon, Yi, & Hartmann 1996). The magnetic stress 
truncates the disk roughly where magnetic pressure halts the slow inward drift of disk material. 
This stress also provides an additional dissipative heating term. Both the truncation radius and 
the heating are essential ingredients for a self-consistent model of the observed disk emission. 
Neither of these quantities can be derived completely from first principles: we thus follow previous 
investigations and parameterize the interaction between the disk and the magnetic field. 

We assume a dipolar magnetic field structure, in which the z-component of the stellar 
magnetic field is 

B(R) « B Z (R) = , (3-1) 

where fi* = B*R% is the magnetic moment of the star, B* is the stellar surface field strength, 
and R* is the stellar radius. The azimuthal component of the induction equation in cylindrical 
coordinates (R, (j), z) is 

dB <t> „ r 9v <t> dv * B <t> , 9 ( n dB A ( o OA 

where rjdiff is the magnetic diffusivity. B$ is amplified by differential rotation in the disk. To 
match the angular velocity of disk material with the star, we adopt a vertical gradient of the 
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rotational velocity, v$ = Rfl, as 



dv ri 



dz 



7 — Q\ 



(3-3) 



where SI* is the stellar angular velocity, Q is the disk (Keplerian) angular velocity, and 7 is a 
parameter which accounts for the uncertainty in the vertical velocity shear (Wang 1987, 1995, Yi 
1995). The amplification dominated by differential rotation is then 



dt 



7 \Q* — Q\ \B Z 



(3-4) 



We model the diffusive loss term by assuming that the magnetic Prandtl number is of order unity: 
Vdiff ~ a (Yi 1995, Campbell 1992). This expression is appropriate when turbulent diffusion is 
responsible for both the diffusive field loss and the viscous angular momentum transport. The 
diffusive loss time scale is Td ~ (ail) -1 and 



dt 



\B(f> 

Td 



aQ \Ba 



In steady state, we equate the amplification and diffusive loss terms: 



B<f) 




7(0, - fi) 


B z 




aft 



(3-5) 
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We set the inner disruption radius, Rq, by balancing magnetic force with the force of disk 
material drifting inward (Wang 1987, Campbell 1992, Yi 1995): 



[BzB^^hR 2 



(3-7) 



The disruption radius lies inside the corotation radius, R c , where disk material corotates with the 
stellar magnetosphere. We evaluate the field components at the disk photosphere, as indicated by 
the subscript "z = H." With this prescription, the disruption radius is computed from 

'PoX 3 / 2 
-Rc, 



1- ^ 



(3-8) 



where R c = (GM^P 2 / '47T 2 ) 1 / 3 is the corotation radius and = 2ir/Sl* is the stellar spin period. 
The dimensionless constant A is (e.g. Yi 1995) 

A = [9.1 x 10 4 ]B2 //j7 i?6 )9M -5/3 p -7/3^-i (3 _ 9) 

where B effJ = ^J^B*/W 7 G, P*, 9 = P*/10 9 cm, M M = M*/M , P^ 2 = P*/WO s, and 
g s" 



M i7 = M/10 17 



-1 



Kenyon et al. (1996) compared the truncation radius derived from equation (3-8) with the 
Ghosh & Lamb (1979ab) model for parameters appropriate to pre-main sequence stars. In general, 
i?o is approximately equal to R c for systems with short rotational periods, low accretion rates, and 
high magnetic field strengths. The disruption radius moves closer to the central star as M and P* 
increase, and as P* decreases. Compared to Ghosh Sz Lamb (1979ab) models, our treatment yields 
larger Rq for high M and high P* systems; we derive smaller Po than Ghosh & Lamb (1979ab) 
when M is low and P* is small. 
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3.2. Magnetic Torque and Disk Heating 



The interaction between the stellar field and the accretion disk produces a torque on both the 
star and the disk. The torque on the star has three parts (Ghosh & Lamb 1979ab): (i) a spin- up 
torque from disk material with R < R c , (ii) a spin-down torque from disk material with R > R c , 
and (iii) an accretion torque carried by gas falling onto the central star. The total torque is 

N - 6 N ° 1-(R /R C )^ (3 " 10) 

where N = M(GM*i? ) 1/2 (e.g., Yi 1995, Campbell 1992; Ghosh & Lamb 1979ab). The 
equilibrium spin N = occurs when Ro/R c = 0.915 (see also Wang 1995). Li (1996) derives 
a smaller value of Rq/R c = 0.7 for equilibrium spin from another treatment of the magnetic 
interaction between the star and the disk. This difference illustrates the uncertainty of 
parameterized magnetic disk models. 

To derive the disk temperature structure, we adopt the vertically averaged versions (Mauche 
et al. 1990) of the momentum equations and energy equation which includes the magnetic torque 
and dissipation (e.g. Shu 1992). The <f> component of the momentum equation is 



where v is the kinematic viscosity coefficient. Mass conservation gives 

M = -2ttRXv r = - / 2irRpV R dz (3-12) 

J—H 

where S is surface density of the disk and H is the half-thickness of the disk. Combining (3-11) 
and (3-12), we get 



M d 
~ 2^R?dR 

Defining the viscous torque 



(m) = ww (" Sfl3 D + h^-"- (3 - 13 » 



N msc = -2ttvZR 3 — (3-14) 



and the magnetic torque 



dNma 9 _ R 2 [B<pBz]z=H = _ R 2l^_Jl B l z ( 3 -l 5 ) 



dR a n 

gives the simple integrated momentum equation 



N msc = M(R 2 Q-R 2 M- [ dR'^L. (3-16) 

JRo d R 
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The total dissipation rate in the disk comes from energy conservation 

where dL/dR is the differential luminosity from an annulus between R and R+dR. For a disk surface 
radiating as a blackbody, the effective disk surface temperature is T e ff = [(dL/dR)/4naR\ 1 ^. 
Combining (3-17) with (3-16), we get 

dR = - Nmsc dR- {n *- n) ^R- (3 ' 18) 

and 

M(R 2 Q - R 2 M = N msc + [ R dR'^L. (3-19) 

To determine the temperature distribution for a specific model, we adopt a set of fixed stellar 
parameters (i?*, M*, P*) and disk parameters (a, 7) and choose the two free parameters (M, P*). 
Observations do not directly constrain the disk parameters. The temperature distribution of the 
magnetized disk is specified by a unique combination a/t/ck-B* rather than individual values of 
P*,a,7 (Kenyon et al. 1996) as long as the disk remains optically thick. We adopt a = 0.3 and 
7 = 1. With the physical variables set, we solve eq. (3-8) for the disk truncation radius, Rq. By 
combining Rq with eqs. (3-15) and (3-19), we derive N v i s and N mag . Finally, the differential 
luminosity distribution dL/dR - and hence the effective disk temperature T e ff oc (dL/dR) 1 / 4 ' - is 
obtained from eq. (3-18). 



3.3. Emission from a Magnetized Disk 



We now apply the magnetized disk model to observations of GK Per. Our model must first 
satisfy the X-ray period constraint, which places the corotation radius at R c = 7.2 x 10 9 cm as 
described in §2. Accretion occurs only if the truncation radius, Rq, lies inside R c , which precludes 
disk models with i?o ^ Rc- The observations indicate that the white dwarf is close to spin 
equilibrium, which favors models with Rq ~ R c as described above. Patterson (1991) estimates 
a mean spin-up rate of 0.0008 sec yr _1 from optical data. The X-ray timing data are consistent 
with this rate, although a constant period can fit the X-ray data alone (Ishida et al. 1992). In our 
model, the white dwarf achieves spin equilibrium when iV = at Ro/R c = 0.915 in eq. (3-10), 
which implies A = 5.856 in eq. (3-8). If we adopt the white dwarf parameters from §2, eq. (3-9) 
requires a simple spin equilibrium curve: 



M17 = 57B 2 effJ 



(3-20) 
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Patterson's (1991) spin-up rate similarly limits the accretion torque to \N\ < 5.9 x 10 34 g cm 2 

s~ 2 for a white dwarf moment of inertia, = 2M (I P 2 /5 = 2.8 x 10 50 g cm 2 . The spin-up line is 
then 

Kl " = (mvw-f B '^ (3 - 21) 

where A is 

A = x 7 l 2 /\l-x^ 2 [ 1 (3-22) 

with x as the solution of the equation 

/2 1-8x^7 _ 064 

l_ x 3/2 ~m 17 - V 26 > 

Figure 2 plots the spin equilibrium and spin-up curves for accretion rates appropriate to GK 
Per in quiescence, M ~ a few x 10 16 g s" 1 (Bianchini & Sabbadin 1983), and at maximum, M ~ 
10 19 g s _1 (§2). The GK Per system should fall between the two curves if our magnetic disk 
model is reasonably accurate. Figure 2 assumes that the accretion disk and disk-magnetosphere 
interaction is in steady-state. The GK Per disk is, however, genuinely time-dependent, and a 
derivation of the white dwarf spin should include the time-variability. A complete theory of 
spin equilibrium for a time-dependent disk is beyond the scope of this paper. We consider the 
time-averaged properties as a compromise. For a recurrence time scale of ~ 1000 d and an outburst 
duration of ~ 50 d, the time-averaged accretion rate is ~ 1 — 3 x 10 18 g s -1 . The observed spin 
equilibrium may correspond to an average accretion rate that is much smaller than the outburst 
accretion rate needed for the UV emission. If we select magnetic parameters appropriate for 
spin equilibrium at this accretion rate, disks in outburst spin up the star and have less magnetic 
heating than disks on the spin equilibrium line (see below). The temperature distribution then 
closely approximates a steady-state disk, which does not explain the UV emission (Figure 1). 

In addition to the coherent X-ray pulsation, GK Per also has optical quasi-periodic oscillations 
(QPOs) with amplitudes of a few per cent. The observed QPO period, Pqpo, depends on state of 
the system. Patterson (1981) first identified Pqpo ~ 380 s during quiescence in 1978; Mazeh et al. 
(1985) later estimated a QPO period Pqpo ~ 360 s during the decline phase from the 1983 data. 
These two estimates are consistent within the observational errors. Mazeh et al. (1985) also found 
Pqpo ~ 400 s on two nights in outburst. 

The QPO periods provide another constraint on the magnetized accretion disk model. This 
constraint is, however, highly model-dependent due to the lack of a unique model for QPOs. In 
the beat frequency model (e.g. Lamb et al. 1985), blobs orbiting at a frequency slightly shorter 
than the stellar rotation period produce QPO at the inner edge of the disk. The Keplerian period 
at R = Rq is 

Po = PqpoP*/{Pqpo + P*). (3-24) 

The observed QPO periods indicate Po ~ 180 s and Rq ps 4.6 x 10 9 cm in quiescence, with Po ~ 
190 s and Rq ~ 4.8 x 10 9 cm in outburst. This behavior is opposite to our expectations, because 
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the inner radius should decrease when M is large (see equations (3-8) and (3-9)). This change 
could indicate slight spin-up, although we cannot make any quantitative conclusions. Nevertheless, 
Fig. 2 plots the curve for Pqpo ~ 400s during outburst. Again using the beat frequency model, 
the — M curve for a given Pqpo is 

M 17 = 3.34 x 10 2 (^Pj Blf f , 7 (3-25) 

where 

which becomes ~ 0.66 for P* = 351s and Pqpo ~ 400s. The relation between the accretion rate 
and the magnetic field strength is then Myj = 36.2B^jj 7 . 

Patterson (1991) considered a QPO model where blobs outside corotation produce periodic 
light variations and estimated Rbiob ~ 9 x 10 10 cm in quiescence and Rbiob ~ 3 x 10 10 cm in 
outburst. Although these large radii are inappropriate for the inner disk radius, their ratio is 
comparable to the radius variation expected in a magnetic disk model when M increases during an 
eruption. We speculate that non-axisymmetric disk features outside the corotation region could 
produce the observed QPOs. X-ray irradiation of local density enhancements in a non-steady 
disk could produce bright optical blobs that give rise to QPOs. These density enhancements 
move radially with time in disk instability models (e.g., Kim et al. 1992) and might produce 
the observed QPO period changes described above. Detailed disk instability models with X-ray 
irradiation would test this hypothesis. 

Having isolated magnetic parameters appropriate for GK Per, we now consider the 
temperature distribution of a magnetically truncated accretion disk derived in a self-consistent 
way from the magnetic torque. This torque has its largest affect on the temperature distribution 
when the disk lies well above the spin equilibrium line and the star is spinning down. The torque 
can modify the temperature distribution significantly on the spin equilibrium line, but the disk 
quickly approaches the non-magnetic temperature distribution when M and -B* act to spin up the 
star. We thus do not consider magnetic heating on the QPO curve in Figure 2. 

Figure 3 shows temperature distributions for 10 models along the spin equilibrium and 
spin- up curves of Figure 2. Each pair of models has the same M and Rq; the spin- up models have 
smaller field strengths and hence less magnetic heating. In all cases, the spin equilibrium model 
has a lower temperature near Rq and a higher temperature outside Rq compared to the spin-up 
model. The temperature gradient outside i?o is fairly independent of the magnetic heating and is 
shallower than T oc R~ 3 / 4 in all of the models. 

Figure 4 plots UV spectra for several magnetic disk models along the spin equilibrium and 
spin-up curves. Although the temperature distributions of these models are flatter than steady 
disks, the model fluxes fit the data poorly. The lowest M models (1 and 2) have nearly the correct 



- 11 - 



UV slope, but the UV fluxes lie a factor of 10-20 below the observations. These spectra are similar 
to the spectra of models DI and SS2 in Figure 1. The predicted fluxes reach the correct luminosity 
for M ~ 10 19 g s _1 (models 7 and 10). These models have shallower UV slopes than steady 
non-magnetic disk models at the same accretion rate, but are too blue compared to observations. 



3.4. Other Magnetized Disk Models 



Our results indicate that magnetic heating models cannot explain the UV spectrum of GK 
Per at maximum light. Before we abandon this model completely, we consider if other treatments 
of a magnetized disk have a better chance of reproducing observations. For instance, Mauche 
et al. (1990) introduced two additional parameters to model the magnetic pitch distribution on 
phenomenological grounds. Force- free models also treat the magnetic stress differently and could 
lead to an "improved" temperature distribution compared to our derivation. 



In the original Ghosh-Lamb model (1979ab), the wind-up of the azimuthal field component is 
given by 

= 7 - fi| |BJ . (3-27) 



dt 

The magnetic field changes on a characteristic reconnection time scale « \vA<j>\ where \vA<j>\ 
is the Alfven velocity associated with and £ < 1 is a numerical reconnection factor. The 
equilibrium field distribution is 

i<>. - <>; 







7 H 


B z 







(3-28) 

where vaz is the Alfven velocity corresponding to B z . Following a procedure similar to §3.2, the 
torque on the star is 



N 
W 



TT = 1 + 



64 (R /R c ) 31 / 16 
21 1 - (R /R c ) 3 / 2 



7 (Ro 
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3/16 
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7 (R, 
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3/16 
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1 (R, 



out 



8 V Rr. 
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The disk truncation radius Rq is set by eq. (3-9) with A replaced by 



A = [4.4 x 10 3 



£ a 9/20 



2/3 p-29/24^-23/40 



*,2 



17 



(3-29) 



(3-30) 



where we have adopted the gas pressure and the Kramers' opacity for the disk structure (e.g. 
Frank et al. 1992). For R^jRc ~ 20, spin equilibrium is reached when Rq/R c = 0.591, which 
indicates that the disk extends further inward from the corotation radius than our present model. 
The difference in the disk temperature distribution becomes most significant outside the corotation 
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radius. The ratio of the magnetic pressure Pb to the gas pressure P g is Ps/Pg ~ (l/0 2 (R/Rc) 2 , 
so magnetic dissipation dominates at R 3> R c - In this region, 

dL „ lJL n ^ K 2^-13/16 (3 . 31) 

and 

T e// oc iT 29 / 64 ~ 22-V2 (3-32) 

which is much flatter than the familiar T e jj oc i?~ 3 / 4 temperature profile. This result is particularly- 
encouraging for GK Per as outlined earlier. We do not further explore this model, however, due 
to the arbitrary prescription for the azimuthal field amplification. The dynamical structure of the 
outer region with Pb > P g is also largely unphysical (Wang 1987). 

An alternative model requires the magnetosphere to satisfy the force-free condition, i.e. 
VB oc B. The magnetic field lines then continuously reconnect to prevent the growth of 
the magnetic stress inside the magnetosphere. Magnetic reconnection, which is the field loss 
mechanism, occurs on a time scale ~ 7 |0* — Q, ss 17* for f2* > Q and s=s Q, for fi* < Q, (Aly & 
Kuijpers 1990). Eq. (3-6) then indicates that the magnetic pitch must lie below a constant of 
order unity. We can thus describe this model by replacing 7/a with a constant jmax ~ 0(1) 
and Q in the denominator of eq. (3-6) with for R > R c where f2* > fl. The resulting torque 
expression produces an equilibrium spin at Rq = 0.950i? c (Wang 1995). For 7 max = 7/0, this 
model is almost indistinguishable from our model in the region R < R c . The disk temperature 
distribution for R > R c is also adequately approximated by our model. 



4. X-ray Heating 



X-ray heating is our final alternative to produce a flat temperature gradient in GK Per's 
accretion disk. In this mechanism, the disk absorbs radiation from the base of the accretion 
column - and perhaps the surrounding white dwarf photosphere - which raises the local blackbody 
temperature. Wade (1988) and Smak (1989) show that this process is unimportant for most CVs, 
because the disks are small and radiation from the boundary layer has a large dilution factor in 
the disk plane. However, the disk in GK Per is several times larger than a typical CV disk, and 
the accretion column more easily illuminates a disk than a boundary layer. 

To establish that this process is viable for GK Per, we consider disk heating near maximum 
in outburst. The accretion rate is then ~ 10 19 g s _1 (e.g., Kim et al. 1992), which produces 
a luminosity of ~ 250 L & at the white dwarf's surface. For simplicity, we assume the central 
star radiates this luminosity uniformly in all directions and parameterize the height of the disk 
photosphere above the midplane as H/R = 0.04i? n (see Warner 1995; Kenyon & Hartmann 1987). 
A truncated disk with parameters appropriate for GK Per intercepts ~ 2% of the accretion 
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luminosity for a standard disk with n = 0.125; the fraction increases to ~ 6% for n = 0.25 and 
~ 19% for n = 0.5 (see Kenyon & Hartmann 1987). If the disk radiates all of this energy in the 
continuum, the extra luminosity from X-ray heating is 5-50 L & . This extra luminosity is close to 
the total UV luminosity at maximum, ~ 20 L (Wu et al. 1989). 

This calculation is overly simplistic, because it does not consider the geometry of the accretion 
column or the details of X-ray absorption in the disk atmosphere. In a standard magnetic accretion 
geometry, gas accreting onto the central star along magnetic field lines should produce two rings 
of emission at latitudes, ±b, measured from the stellar equator (see Patterson 1994 and references 
therein). Standing shocks form above each ring at a height (e.g. Yi et al. 1992) 

H s = [2 x 10 8 cm] So Mf/X 7 4/ XrMr 7 5/7 (4-1) 

where B*j = f?*/10 7 G and s is a constant of order unity. This approximation is valid even when 
the magnetic axis is not aligned with the rotational axis, although the rings are then only partially 
filled with emission (Konigl 1991; Hellier 1993). This height is small compared to the stellar 
radius at high accretion rates, M > 10~ 18 g s -1 , and we can reasonably approximate the accretion 
hot spots as a point displaced from the center of the white dwarf. Our simple estimates for the 
accretion luminosity intercepted by the disk are then valid. However, the height of the accretion 
shock increases as M decreases and becomes H s ~ i?* for M < 10 18 g s _1 . We then cannot simply 
approximate accretion hot spots as a point close to the white dwarf photosphere, which greatly 
complicates our calculation of disk heating. 

A direct calculation of X-ray processing by the disk is also complicated. Previous results 
suggest that the disk is very optically thick to hard X-rays; this radiation heats the disk to 
produce extra continuum and line emission (Raymond 1993 and references therein). Raymond's 
(1993) models for Sco X-l and other luminous X-ray binaries show that X-rays penetrate close 
to the disk photosphere at small disk radii, R ~ 10 8 to 10 9 cm, where electron scattering is the 
dominant opacity source. These disk radii emit approximately as blackbodies. Line emission is 
more important at large disk radii, R ~ 10 10 to 10 11 cm, where the X-ray optical depth exceeds 
unity well above the disk photosphere. Prominent H I Balmer and He I, He II emission continua 
also modify the blackbody spectrum of the outer disk in Raymond's (1993) models. 

The X-rays in GK Per are a factor of 10 or more weaker than expected from the UV 
luminosity, so we might expect considerable EUV heating of the disk. Proga et al. (1996) 
considered a red giant atmosphere illuminated by a hot white dwarf with T e ff = 2 x 10 5 K and 
showed that the atmosphere re-emits most of the incident radiation as emission lines and a strong 
recombination continuum. The illuminated atmosphere emits 25%-50% of the incident flux as 
continuum radiation for conditions less extreme than those in GK Per. We thus expect a large 
contribution from H I Balmer continuum radiation if the accretion column emits mostly in the 
EUV. This radiation would flatten the UV continuum considerably compared to optically thick 
disk models (see below). 

To place better limits on the importance of X-ray heating in GK Per, we assume a point-like 
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X-ray source illuminating a disk whose photospheric height varies as R n . The X-ray heating flux 
at a disk radius R is roughly 

where 77 < 1 is the uncertainty of the X-ray efficiency at the central source (e.g. Frank et al. 
1992) and / is the X-ray absorption efficiency that takes into account the X-ray albedo ~ 0.3 
and the X-ray heating geometry. The factor / also includes the radial variation of the projected 
surface area of the disk and thus slowly varies as a function of R (see below). This heating source 
dominates locally generated accretion energy for disk radii, 

R>R X = [1.5 x 10 n cm]i?* i9 (^2) ■ ( 4 ~ 3 ) 

Thus heating begins to dominate accretion near the outer edge of the disk in GK Per, 
Rout ~ 1-5 x 10 11 cm, for frj ~ 0.1. The temperature profile is then intermediate between the 
T oc i?~ 3 / 4 of a steady disk and the T oc Rr^l 2 limit of a pure reprocessing disk (see Smak 1989, 
Kenyon & Hartmann 1987, and references therein). For simplicity, we model F x as a local heating 
term and compute the radial temperature gradient in the blackbody approximation: 



T eff = 



1/4 

(4-4) 



Several factors complicate the direct calculation of / needed to derive a self-consistent 
temperature distribution for an X-ray heated disk. In addition to the geometry of the accretion 
column, the white dwarf shields some portions of the disk from hard X-rays when the magnetic 
and rotational axes are not aligned. If the angle between these two axes is (ft, the innermost disk 
radius exposed to X-rays is approximately 

R i = R * (^~J ~ ^ ■ ( 4 ~ 5 ) 
\sm<f) J 

In GK Per, the region near the truncation radius at R ~ Rj ~ 10-R* is shielded for <p ~ 6°, which 
is comparable to the misalignment angles measured in other CVs (see Hellier 1993 and references 
therein). 

We compute / for R > Rj assuming the disk photosphere absorbs and reradiates incident 
X-rays at T e ff from equation (4-11). We adopt an X-ray albedo of 0.3 and the Shakura-Sunyaev 
variation of disk thickness with radius, H/R (Frank et al. 1992; Warner 1995) and estimate 

/ « 2.1 x ID"* (^) M^X/ 20 < (4-6) 

for H/R 1, where R±o = i?/10 10 cm. For GK Per specifically, 

/ w 2.2 x lO^M?/ 20 ^ 8 . (4-7) 
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We note that the dependence on a in eq. (4-6) is very weak. Given the uncertainties, we 
parametrize the X-ray heating as 

/ = hoRio (4-8) 

where f w denotes / = 2 x lO -2 ^ 20 at R = 10 10 cm. This expression recovers equation (4-7) for 
n = 1/8. 

Previous studies of disk heating show that 'flat' disks with n < 1/8 have temperature 
distributions similar to a steady disk with T oc i?~ 3 / 4 (Kenyon Sz Hartmann 1987). Figure 5 shows 
temperature distributions for four combinations of (M, B e ff) that lie between the spin equilibrium 
line and the spin up line at each M in Fig. 2. In each panel, the flattest temperature distribution 
corresponds to n = 1.2, the middle one to n = 0.6, and the steepest one to n = 0.3. 

Figure 6 compares spectra of disks with a low X-ray heating rate (n = 0.3) and disks without 
X-ray heating. The extra luminosity of the heated disks is considerable, even when the disk 
height increases slowly with radius as in standard hydrostatic disk models. However, models 
with the correct UV slope are underluminous compared to observations. The UV slopes of more 
luminous disks are too blue compared to the data. The agreement between the models and the 
data increases with increasing n. Figure 7 compares spectra for two other model temperature 
distributions - model 27 and model 28 - with observations of GK Per. Both models match the 
data over 1750-3250 A and are a factor of ~ 2 above the data at shorter wavelengths. Both models 
predict V magnitudes within ±0.3 mag of those observed near maximum. Less extreme heating 
models also predict the observed V but fail to agree with the UV data over as large a wavelength 
range. 

Although we find some agreement between models and observations with X-ray heating, the 
resulting variation of disk height with radius is severe. A disk with n = 1.2 is extremely concave 
and physically unrealistic. Previous calculations of disk heating indicate that the disk probably 
maintains hydrostatic equilibrium with n < 0.3-0.5 (e.g., Kenyon & Hartmann 1987; Raymond 
1993 and references therein) and cannot reach n > 0.5. Even if an n = 1.2 disk can be achieved in 
nature, the outer rim of this disk occults most of the inner disk for modest inclinations, i « 45°. 
Thus, this disk model is probably inappropriate for GK Per. 

Finally, we consider a simple optically thin disk model. We adopt the X-ray heating model 
described above for n = 0.3, assume X-ray photons ionize material above the disk photosphere, 
and calculate the recombination continuum for an LTE slab of hydrogen with an optical depth 
at 3000 A, T3000 ~ 0.5. We adopt a slab temperature of 15,000 K, which results in a total UV 
luminosity, 5 L , comparable to the accretion energy intercepted by the disk. Although these 
parameters are not 'best' for an accretion disk corona in GK Per, this calculation provides a 
first approximation to optically thin emission in the heating model. Figure 8 compares spectra 
of spin equilibrium model 8 with and without a recombination continuum from our accretion 
disk corona. The recombination continuum flattens the UV spectrum considerably, although the 
predicted spectrum falls below the observations at all UV wavelengths. This result is sensitive to 
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our adopted parameters. Hotter coronae produce very blue spectra; smaller optical depths emit 
too little UV flux. A better calculation, similar to Raymond's (1993) treatment, is needed to show 
if this hypothesis can explain the spectrum in detail. 



5. Summary and Discussion 

We have investigated several accretion disk models to account for the observed UV spectrum 
of GK Per. This spectrum indicates a large inner disk radius and an unusually flat temperature 
distribution. These requirements are not satisfied by any non-magnetic steady disk model or by 
existing disk instability models that are otherwise successful in explaining the recurrence time 
scales of GK Per's occasional dwarf nova outbursts. 

Magnetic disk models also fail to explain the UV spectrum. We used data for QPOs and 
the spin-down rate to constrain the inner radius of an equilibrium system and found that M and 
-B* are confined to a narrow band with 36 < Mn/B* ff 7 < 57. Magnetic torque significantly 
affects the disk temperature distribution at the upper end of this range and produces shallower 
temperature laws than predicted for standard steady disks. Our model temperature distributions 
are, however, steeper than needed to explain the observed UV spectrum. 

X-ray heating models are a promising mechanism to produce a flat disk temperature gradient 
in GK Per. Optically thick models that absorb and reprocess X-rays close to the disk photosphere 
require significant curvature in the disk to reproduce the UV spectrum. These disks are physically 
unrealistic. Optically thin models can explain the UV spectrum if Balmer continuum radiation 
dominates the emitted spectrum of an accretion disk corona. Our example in Fig. 8 represents a 
favorable case where the disk+corona intercepts 5%— 10% of the accretion luminosity and remains 
cool enough to emit much of this energy at 2000-3000 A. More detailed calculations are needed to 
test the viability of this solution. 

The heating model has several consequences that can be directly tested by observation. If 
much of the 2000-3000 A continuum is H I Balmer emission, this flux should show large amplitude 
variations phased with the X-ray light curve. Patterson (1991) detected stable variations at U 
with the appropriate period. The amplitude of this variation is small, ~ 0.004 mag, compared 
to the X-ray amplitude. Observations of emission shortward of the Balmer jump would provide 
a better test of the heating model. The optical emission line spectrum provides another test. 
Measurements of the H I line fluxes through an outburst directly constrain the flux from the 
Balmer continuum; detection of periodic variations in line flux phased with the X-ray period 
places additional limits on the heating model. 

Even if X-ray heating does not explain the UV spectrum directly, this heating is an important 
contribution to the energy balance of the disk itself. Our estimates in §4 indicate that the disk 
intercepts ~ 5% of the accretion energy radiated by material impacting the surface of the white 
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dwarf. This energy is a significant fraction, ~ 10%-20%, of the disk luminosity in outburst and 
probably changes the vertical temperature structure of the disk. This heating could be larger 
during quiescence when the height of the accretion shock above the white dwarf increases (cf. eq. 
(4-1)). The vertical temperature structure is an important feature of disk instability models that 
trigger outbursts when disk material is too cold to radiate viscous energy generated locally (see 
Kim et al. 1992 and references therein). Thus, X-ray heating might modify the properties of dwarf 
nova-like behavior in this system. 

The accretion rates required for the observed UV emission at maximum exceed 10 19 g s" 1 and 
are incompatible with the observed X-ray luminosity. The observed X-ray luminosity (L x ) implies 
M ~ L X R*/GM* < W 18 g/s (Ishida et al. 1992, Norton et al. 1988, Yi et al. 1992, Yi 1994, Yi & 
Vishniac 1994) in the X-ray emitting magnetic accretion column. The X-ray luminosity remains 
more than an order of magnitude higher than the optical/UV luminosity in quiescence, which 
further contradicts the results of disk instability calculations (Yi et al. 1992, Kim et al. 1992). 

In a magnetically channeled accretion column (Frank et al. 1992), the radially free- falling 
material shocks above the white dwarf photosphere. Highly ionized material in the postshock gas 
emits X-ray through bremsstrahlung. Due to highly anisotropic (Compton) scattering geometry, 
most photons travel radially outward along the column axis and exert a radiation drag on the 
infalling material. As a result, the X-ray emission temperature can be substantially lower than the 
expected postshock temperature ~ 5 x 10 S K. This effect becomes significant for high M and small 
column cross-sectional area (Yi 1994, Yi &: Vishniac 1994). For typical parameters derived from 
the UV spectrum and spin, the ratio of the column cross-sectional area to white dwarf surface 
area is ~ 0.06s o (M/10 18 5/s) 2 / 7 , where s Q < 1 is a constant (Yi 1994). The radiation-modified 
postshock (X-ray emission) temperature falls below 10 8 K when M ~ 4 x 10 18 g s _1 for s Q ~ 0.5 
(Yi 1994, Yi& Vishniac 1994). The X-ray temperature continues to fall as M increases; the 
X-ray luminosity is then much less than the total accretion luminosity. For example, the X-ray 
luminosity for E > 2 keV at M ~ 4 x 10 18 g/s implies an apparent M ~ 6 x 10 17 g s" 1 . The 
predicted X-ray emission spectra (after photoelectric absorption) tend to show little indication of 
this effect (Yi & Vishniac 1994). We conclude that the mass accretion rates in UV and optical 
can be reconciled with the observed X-ray luminosity. This effect requires M > 10 18 g/s which is 
consistent with the observed UV emission. 

Finally, combining spectral information and improved spin torque measurements could 
convincingly constrain physical parameters which are largely unknown. Precise measurements 
of spin-up (-down) torque during the rise (decline) phase of outburst could provide especially 
valuable information on the magnetized accretion disk model. More reliable data on QPO during 
quiescence and outburst could further constrain the magnetized accretion disk model. 

We thank D. Proga, J. Raymond, and J. Craig Wheeler for helpful advice and discussions. I. 
Y. acknowledges support from SUAM Foundation. S. K. acknowledges support from NASA grant 
NAG5-1709. 
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Fig. 1. — Ultraviolet spectra of the two steady disk models and the Kim et al. (1992) disk instability 
model (see text). Filled circles plot dereddened IUE fluxes for spectra acquired at maximum during 
an outburst (SWP 13497 and LWR 10143 on 1981 March 15; Wu et al. 1989). For inclinations 
i > 45°, these models cannot fit the spectral slope and the flux level simultaneously. 
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Fig. 2.— Magnetic field (B effJ = v / ^JaB i ,/10 7 G) - mass accretion rate (M 17 = M/W 17 g/s) 
parameter space. The long dashed line corresponds to the observed QPO period of Pqpo = 400s 
in the beat frequency model. The short dashed line corresponds to the spin-up torque derived in 
Ishida et al. (1992). The dotted line corresponds to spin equilibrium. The parameter space above 
the spin equilibrium line - the spin-down region - is ruled out from the observed spin periods. Ten 
models marked by numbers are considered for the standard magnetized accretion disk model. The 
temperature distributions corresponding to these models are shown in Figure 3. Models numbered 
11-28 are X-ray heated models shown in various other figures. 



-22- 




9.5 10 10.5 11 11.5 9.5 10 10.5 11 11.5 




9.5 10 10.5 11 11.5 9.5 10 10.5 11 11.5 

log R(cm) 




9.5 10 10.5 11 11.5 
log R(cm) 



Fig. 3. — The temperature distributions for the ten models in Figure 1. The disk temperature 
distribution is set by the mass accretion rate and the magnetic torque heating. In each panel, the 
solid line corresponds to spin-up and the dotted line to spin equilibrium. Higher magnetic fields 
for spin equilibrium raise disk temperature through magnetic torque heating and shifts the disk 
temperature maximum to larger radii. 
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UV spectra of some of the spin-equilibrium and spin- up models shown in Figure 2. 
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Fig. 5. — Disk temperature distributions for X-ray heating models. The models are all in slight 
spin-up allowed by the existing data. The parameters are shown on each panel. In each panel, 
three different X-ray heating models are shown (solid line: n = 1.2, dotted line: n = 0.6, dashed 
line: n = 0.3). The locations of the models in Figure 2 are marked by the model numbers 11-22. 
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Fig. 6. — Comparison of spectra of disks with low X-ray heating (n = 0.3) and disks without X-ray 
heating. The model numbers are identical as those in Figure 2. 
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Fig. 7. — UV spectra of two high X-ray heating models 27 (n = 1.3) and 28 (n = 1.4). These good 
fits are due to the very flat temperature distributions resulting from a large X-ray heating rate as 
described in the text. 
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Fig. 8. — Comparison of spectra of the spin equilibrium model 8 (n = 0.3) with and without 
additional continuum from a simple accretion disk corona. Models with a recombination continuum 
fit the data but are uncertain as described in the text. 



